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Theoretical results on the dynamics of dislocations in Rayleigh-Benard convection are reported 
both for Swift-Hohenberg models and the Boussinesq equations. For intermediate Prandtl numbers 
the motion of dislocations is found to be driven by the superposition of two independent contribu- 
tions: (i) the Peach-Koehler force derived from the change of a Lyapunov potential with pattern 
wave number; (ii) a non-potential advection force on the dislocation core by its self-generated mean 
flow. Their competition allows for the first time to understand the experimentally observed bound 
dislocation pairs. 

PACS numbers: PACS: 47.54.+r, 47.20.Ky, 47.20.Lz, 



Striped patterns are ubiquitous in nature. They 
are found in physical, chemical, and biological systems, 
which are driven away from equilibrium jjj. In general 
natural patterns are not perfect due to the presence of 
defects, like grain boundaries (line defects) and topolog- 
ical point defects (dislocations) 0- The nucleation, mo- 
tion, and annihilation of dislocations is essential for many 
pattern-selection processes, which are initiated by modu- 
lational instabilities Q . Dislocations govern the ordering 
kinetics of initially disordered patterns (^] and sustain 
in defect turbulent systems the perpetual reordering of 
the planforms Q. Thus much effort has been devoted 
to the study of dislocations in striped patterns (see e.g. 
0) §|> [| 01)- Dislocations present a simple realization of 
topological singularities in a field description of continu- 
ous extended systems M . To which extent their dynamics 
can be understood in terms of "particles" subject to effec- 
tive forces is a general important issue, that transcends 
pattern formation. 

Here we study dislocations in Rayleigh-Benard con- 
vection (RBC), which in the past has proven itself as a 
paradigm for pattern forming systems ||, ||. In RBC 
a horizontal layer of a simple fluid is heated from be- 
low and cooled from above. For temperature differences 
AT above a critical value AT C , buoyancy driven convec- 
tion sets in typically in form of striped convection roll 
patterns with the critical wavenumber q = q c . The ba- 
sic hydrodynamic equations (Oberbeck-Boussinesq equa- 
tions, OBE) are well established and the theoretical anal- 
ysis can be quantitatively compared with well controlled 
experiments The communality of RBC with other 
bulk pattern-forming systems like, for example, vibrat- 



ing granular layers gas discharges jjjj, and Taylor- 
Couette flow |jj is most clearly expressed in universal 
model equations, like the Swift-Hohenberg (SH) equation 



jl2[ and its generalizations 13 



In this letter we concentrate on the dislocation dy- 
namics in the case of gas convection, which over the past 
decade has attracted most interest |J. Besides the ap- 
propriate SH-model, for the first time the full OBE equa- 
tions are analyzed. We found the motion of dislocations 
to be driven by two independent "forces". One is the 
well known Peach-Koehler force, which describes the ten- 
dency of the system to develop towards a striped pattern 
with an optimal average wavenumber q « q c || |(|. Of 
particular importance is our explicit identification of a 
second non-potential force due to the advection of the dis- 
location core by its self-generated, roll-curvature driven 
mean flow. We demonstrate that this force acts in gen- 
eral to remove dislocations from the system as to reduce 
q. The two forces may balance each other at a certain 
background wave number q = qo < q c , such that the dis- 
location becomes stationary. For patterns with q > qrj 
the competition of the two forces can prevent the an- 
nihilation of two dislocations with opposite topological 
charge. This explains for the first time the experimen- 
tally observed bound dislocation pairs |H . Although the 
significance of the mean flow had been expected for a long 
time in line with experiments ||l5f , its direct identification 
in a previous theoretical analysis failed |l6| . 

In RBC the strength of the mean flow is determined 
by the Rayleigh number R (a dimensionless measure of 
the temperature gradient across the convection cell) and 
in particular by the Prandtl number a = v j k, with the 
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kinematic viscosity v and the thermal diffusivity k. In- 
termediate Prandtl numbers (p ~ 1) are realized in gas 
convection experiments In Fig. [I], a snapshot of the 
midplane temperature field tp(x,y) containing a disloca- 
tion with positive topological charge Q = 2tt from simu- 
lations of the OBE is shown. Q is defined via the winding 
number of the phase gradient V<I>(a;, y) oitp around a dis- 
location with the two possibilities § VQ(x,y)ds = ±2ir. 
The mean flow (arrows), which is driven by strong cur- 
vature of the rolls in the vicinity of the dislocation core, 
is maximal along the symmetry axis and tends to advect 
the dislocation downwards out of the system. In the di- 
rection perpendicular to the plane shown in Fig. ^ the 
mean flow shows an almost perfect parabolic profile. 




FIG. 1: Midplane temperature field (gray scale) of a disloca- 
tion with positive topological charge Q = 2ir and the super- 
imposed mean-flow velocity field (arrows). The dislocation is 
moving downwards along the y -axis. From simulations of the 
OBE for e = (AT — AT C ) / AT C = 0.3, a = 1.2 and background 
wavenumber q = q c (see text). A dislocation with Q — — 2n 
would move upwards in this case. 

While a comprehensive exploration of the i?, a param- 
eter space within the OBE is extremely time consum- 
ing, we have gained much insight into the basic physi- 
cal mechanism of dislocation dynamics by first studying 
the standard two-dimensional generalized SH-cquations 
(GSH) 
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In terms of the two-dimensional field i/j(x,t) the GSH- 
equations describe (in suitable dimensionless units) the 
bifurcation to patterns with critical wavenumber q c = 1 
and their nonlinear saturation. £5 oc e = (AT/AT C — 1) 
measures the relative distance to onset. The coupling 
to the mean drift (or mean flow) U(x) is characterized 



by a coupling constant g (in RBC g cx 1/er) and a cut- 
off parameter c. For g = the GSH-equations derive 
from a Lyapunov functional C JIJ and any dynamics is 
"downhill" towards the minimum of C. 

In our numerical analysis we focused on dislocation 
climb, i.e., the motion parallel to the roll axis as shown in 
Fig. 1, in a rectangular domain —L XiV /2 < x, y < L XjV /2 
with width L x and length L y = 2L X . We solved both the 
GSH- and the OBE equations numerically by a pseudo- 
spectral method with semi-implicit time stepping [fP7| . 
In order to minimize finite-size effects we have simu- 
lated dislocation pairs with periodic boundary condi- 
tions in the x, y— directions (Fig. 1 has to be extended 
mirror-symmetrically along the y-axis). We found, how- 
ever, almost identical results for the less expensive sim- 
ulations with a single dislocations in a box kept finite 
in the y direction by gradually ramping e to zero at 
y = ±L y /2. The ideal pattern without an immersed 
dislocation consisted of up to N — 64 roll pairs of a 
wavelength A = 27r/q_, corresponding to L x = NX. The 
typical numerical resolution was > 8 grid points per roll 
diameter. It is crucial to study systematically the de- 
fect dynamics as function of the background wavenum- 
ber q of the underlying pattern. A superimposed dis- 
location leads to N + 1 roll pairs in some parts of the 
cell (see Fig.l, below the dislocation) and consequently 
to a reduced wavelength A + = N/(N + 1) A and thus 
to a wavenumber q + — 2ir/X + > g_. We found that 
defining the background wavenumber as q = (q + + q_)/2 
was most effective to absorb the finite-^ corrections. In 
designing the details of our calculational scheme and to 
validate our SH-code we have considerably profitted from 
the comparison with Ref . || , devoted to the SH-model. 

First we discuss the solutions of the GSH-equations. 
The simulations (for simplicity we consider the case of 
a single dislocation which climbs along the y-axis) were 
initialized with an approximate ansatz tf)(x,y,t — 0) oc 
cos[q- x — 4>(\r\)] at the center fo = (xo,yo) — (0,0) of the 
simulation domain. Here <j>(r) is the polar angle about tq. 
After the initial transients died out (say at \y ~yo\ > 4 A 
in Fig.|l|) the dislocation would climb with a constant 
velocity v st and a relaxed asymptotic shape tp(x,y,t) = 
ip(x,y — v s tt) was reached. In subsequent runs, starting 
from this state considerably reduced the computational 
time. In agreement with prior investigations ^ we 
have identified a wavenumber gu(e, g), such that for q > 
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qo the dislocation climbed downwards (the case shown in 
Fig. 1) with v st < while it moved upward for q < qo- 

The key for the physical interpretation of the results 
lies in a certain universality of ip and U. As demonstrated 
in Fig. |^A for a representative example of uniformly 
climbing dislocations at fixed background wave number 
q = 0.98 the numerical solutions i/j(xo,y) = ipo(y) with 
y = (y — v s t t) of Eq. ^ are virtually identical along the 
symmetry axis (x — xq — 0) although the mean-field 
coupling g varies considerably (2 < g < 20). In addition, 
as shown in Fig. ||B for the same range of parameters, 
the (suitably reduced) mean flow component U y (xo,y)/g 
along the symmetry axis has practically a ^-independent 
shape Uo(y), while the transverse component U x vanishes 
for symmetry reasons. Inspection of Eq.([j]) suggests that 
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FIG. 2: Profiles of U Q {y) = U y (x ,y)/g and ^o{y) = ip(xo,y) 
plotted along the roll axis centered about the dislocation core 
at (xo = 0,y = y — v s t — 0) obtained from SH-simulations 
at es = 0.3 and q = 0.98 in the range 2 < g < 20 for N = 
16, q~ — 0.95 and c = 2. The profiles show approximately a 
fifth of the calculational grid of 512 pixels along y, where 16 
grid points correspond to length of 3 A = 6iv/q-. 

ip(xo, y) can be well approximated by the potential solu- 
tion i/>o(y) a t 9 = 0. This applies also to the source term 
in Eq. (||) with respect to its ^-dependence. Thus Eq.[j] 
on the symmetry axis reduces to: 
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Averaging Eq.([|) over y leads to the following relation 
for v s t as function of q and g: 



v st {q,g) = v pot {q,g^0) + c 1 {q)gU o 
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With respect to the advection term, we have made ex- 
plicit the (reduced) mean flow Uq 1 at the dislocation core 
(i.e. the minimum of Uo(y — 0) in Fig. |^B), which dom- 
inates the average in Eq. ([|) also because of the strong 
maximum of dy i/>o (y) at y — (see Fig. A) Q . The 



quantities v po t , ci and Uq have to be determined numer- 
ically as function of q and es- 

For clarity the consequences of Eq. (||) are discussed 
for Q = 2tt (see Fig. [I]), since the case Q = —2ir is 
analogous. The potential or Peach-Kochler contribution 
Vpot , originating from the functional derivative of C in Eq. 
(^), is negative for q > q c = 1, vanishes for q w q Cl and 
becomes positive for q < q c . Thus the dislocation climbs 
as to move the wavenumber of the pattern in its wake 
towards an optimal value q op t ~ q c , where the potential 
energy has its minimum || . Since c\{q) in Eq.(^|) was 
found to be always positive for all q and es, the mean 
flow (oc Uq 4 < 0) would always advect the dislocation 
such as to decrease the wavenumber q. In fact, for g ^ 
in patterns with q in the range qo < q < q pt the 
potential contribution v po t is overcome by the mean flow 
contribution leading to v s t < 0. At q — qo(e,g) both 
effects cancel leading to a immobile dislocation (v st — 0) . 
In agreement with prior investigations [^| our numerical 
results show that v s t varies almost linearly as function 
q according to v s t oc (q — qjj) over a wide q-range for 
fixed g > and e s > 0. According to Eq.(^) this implies 
at fixed q a linear dependence of v s t on 9, which was 
numerically confirmed in a wide <?-range (0 < g < 30). 

Although in full hydrodynamics (OBE) a direct analog 
of potential C in the GSH is not defined, the reasoning 
in terms of a potential force ( a — + oo) and a competing 
mean flow force seems to apply well. As a typical example 
for a = 1.4, we show in Fig. ||A the linear variation of 
v st (q) as function of q, which crosses v s t = at q = 
qo = 2.75. In Fig. ||B it is demonstrated, that v s t varies 
like the SH-results linearly in o~ l ~ g at fixed q. Our 
findings are in perfect agreement with the experimental 
and numerical analysis of the dynamics of the off center 
" giant spirals" for crwl where the uniformly rotating 
outer tip would probe roll patches in a fairly large q— 
range. 

In agreement with experiment |l4j , from equation (|5|) 
the possibility of bound dislocations pairs is obvious. A 
well separated pair of oppositely charged dislocations on 
a background pattern with qu < q < q Cl is driven by 
the dominant mean-flow contribution oc |£7g^| in Eq. (|^) 
towards their annihilation. However, this advection force 
weakens continuously, as the opposing mean flow contri- 
butions from each dislocation begin to overlap destruc- 
tively. Eventually the dislocations stop moving, when the 
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FIG. 3: Dislocation velocity v st in units of 10~ 2 d/t v with d 
the cell width, and t v = d 2 / 'k the vertical diffusion time. A: 
as function of q at fixed a — 1.4. B: as function of a in the 
range 0.8 < a < 2.5 of at q — 2.91, data from simulations of 
the OBE for AT/AT C - 1 = 0.3. 

reduced advection force is balanced by the opposing po- 
tential force. To our surprise even the deceleration of the 
approaching dislocations was well captured adiabatically 
by Eq. (^]), derived under the assumption of constant v s t- 
In fact, we found a direct proportionality between the de- 
creasing velocity v(t) of the approaching dislocations and 
the resulting Uq* (t) . The equilibrium distance I d of the 
dislocations in the bound pair as function of q and <?, or 
<7, respectively, has not been investigated systematically. 
However, typically we find Id ~ 2 — 3A Jl9[ . 

In our OBE and GSH-simulations dislocations are 
found to be very robust objects except when q lies in 
the vicinity of a cross-roll (CR) instability boundary |2l[] 
for ideal rolls, e.g. in the OBE near q — 2.48 < q c = 3.11 
for g = 1.4, e = 0.3. Then, as in the experiments H |2(| 
the dislocation core might split into a chain of bubbles 
by a localized transverse bridging of adjacent rolls. This 
local CR instability has been described at first in the 
SH-model simulations p2[ , but we doubt that the discli- 
nation mechanism invoked there has a merit. 

In conclusion, we have presented an analysis of the de- 
fect dynamics for medium a and finite R in RBC, which 
has allowed to separate clearly potential- and mean-flow 
forces. Our main interest was to elucidate the basic 
mechanism. Details, like for instance the dependence of 
qjj on o and e and the impact of the stability proper- 
ties of the background pattern will given in a separate 
manuscript. We expect that our methodology, to sepa- 
rate potential and nonpotential effects will be applicable 
as well to a number of other systems that show striped 
patterns as mentioned before. E.B. is grateful for sup- 
port from the National Science Foundation under grant 
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